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ABSTRACT 



We construct a holographic map from the loop equation of large- N QCD 
in d = 2 and d = 4, for planar self-avoiding loops, to the critical equation 
of an equivalent effective action. The holographic map is based on two in- 
gredients: an already proposed holographic form of the loop equation, such 
that the quantum contribution is reduced to the evaluation of a regularized 
residue; a new conformal map from the region encircled by the based loop to 
a cuspidal fundamental domain in the upper half-plane, such that the regu- 
larized residue vanishes at the cusp which is the image of the base point of 
the loop. The critical equation of the holographic effective action determines 
a unitary Abelian local system in d — 2 and a non-Abelian twisted local 
system in d = 4. As a check in the d = 2 theory, we study the distribution of 
eigenvalues of the Wilson loop implied by the critical equation. As a check 
in the d = 4 theory, we study the first coefficient of the beta function implied 
by the holographic loop equation and, as a preliminary step, that part of the 



second coefficient which arises from the rescaling anomaly, in passing from 
the Wilsonian to the canonically normalised (holographic) effective action. 



1 Introduction 



This paper grew out of the attempt to obtain directly from the loop equation 
[1, 2] for self-avoiding loops the celebrated formula for the distribution of the 
eigenvalues of a Wilson loop of area A in the weak coupling phase, in the 
large- N limit of two dimensional QCD on a sphere, for g 2 A small with respect 

to g Agpfoer,,'. 

r N 

w(A)= /n^ni^-^iexp(-^£^ 2 )x 

x^N^coaBi (1) 

i 

An essentially equivalent formula was first obtained from the Wilsonian lat- 
tice action using functional techniques [3] and afterward Eq.(l) was derived 
from the "heat kernel action" [4, 5] or from the semi-circle law for the eigen- 
values of free random variables [6] , which are employed in the complete oper- 
atorial solution of the loop equation in d — 2 for arbitrary self-intersections 
[7]. We present here a derivation based on the loop equation restricted to 
self-avoiding loops but that extends to the (planar) four dimensional case. 
Restricting to (planar) self-avoiding loops has some virtue, since the general 
solution for arbitrary self-intersection seems to be out of the reach of our 
methods in d — 4 and in d = 2 as well. In this respect, already long ago, 
the loop equation in d = 2 and in d = 4 has been written in terms of the 
distribution of the eigenvalues for self-avoiding loops [8] 1 . In Eq.(l) and 
in this paper we have ignored the contribution of the exterior of the loop, 
assuming that the internal region is much smaller than its complement. In 
addition we have not taken into account the periodicity of the eigenvalues. 
We believe that the techniques of this paper can be extended to include also 
the more general case. However for simplicity we will consider it elsewhere. 
In our derivation there are two basic ingredients. The first one is a holo- 
graphic form of the loop equation for planar loops in d = 2 and d = 4. This 
holographic loop equation was derived in [9]. Using by analogy the language 
of the correspondence between the boundary Af = 4 four dimensional gauge 
theory and the bulk five dimensional super-gravity (string theory) [10], by 

1 We thank the referee for pointing out ref. [8] to our attention. 
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holography in this context we mean a correspondence between the loop equa- 
tion, that we think as a theory defined on boundary curves, and an equivalent 
effective action for the eigenvalues of a Wilson loop, that we think as a the- 
ory defined on the bulk. The holographic form of the loop equation that 
we refer to is a preferred form of the loop equation obtained by means of 
appropriate changes of variable so that this boundary-bulk correspondence 
becomes almost manifest, for the reason that the quantum term in the loop 
equation, given as usual by a contour integral along the loop, is computed as 
a regularized residue, loop independent for self-avoiding loops. This loop in- 
dependence makes the holographic loop equation close to admit an equivalent 
(holographic) effective action. To this form of the loop equation we associate 
a "classical" holographic effective action, r. T is "classical" in the sense that 
the loop equation for T still contains a non-vanishing quantum contribution 
given by the regularized residue. Yet, T already contains quantum correc- 
tions and indeed it is related to the quantum holographic effective action, 
T q , by means of a suitable conformal mapping and gauge fixing. In fact the 
second ingredient needed to complete the holographic correspondence is a 
conformal map of the region encircled by the based loop that occurs in the 
loop equation to a cuspidal fundamental domain. On such domain an ef- 
fective action for the distribution of the eigenvalues of the Wilson loop, T q , 
can be indeed constructed, because the regularized residue vanishes. Hence 
the holographic loop equation becomes equivalent to a critical equation for 
T q . The second ingredient, i.e. the conformal map, is essentially new in this 
paper with respect to [9] although it already appeared in implicit form in 
appendix C in [9]. Yet, there, the two-steps logic of holography and confor- 
mally mapping was somehow mixed so that we did not realize that the two 
dimensional technique of appendix C could be in fact extended to the four 
dimensional case as well. More precisely, our quantum holographic effective 
action furnishes a critical equation for the eigenvalues of the curvature of 
a unitary Abelian local system [11, 12, 13, 14, 15] in the two dimensional 
case and for the eigenvalues of the curvature of a non-unitary non- Abelian 
twisted local system in the four dimensional case. By a twisted local system 
we mean here a central extension of a (possibly infinite dimensional) repre- 
sentation of the fundamental group of a punctured Riemann surface. Thus, 
while the curvature of the local system determines directly the eigenvalues 
of the Wilson loop in d — 2, it does it only indirectly, via the non- Abelian 
gauge connection, in d = 4. 
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To summarise, the basic ingredients of our holographic correspondence are 
the following ones. 

A holographic form of the loop equation in d — 2 and d = 4 as described 
in [9]. This form of the loop equation follows from changes of variable, that 
though have a geometric meaning in terms of symplectic reduction to a mi- 
crocanonical ensemble, followed by a choice of a holomorphic gauge, do not 
involve in a crucial way algebro-geometric concepts. These are rather needed 
in the assignment of a local system with a lattice of punctures in the d — 2 
theory and of an infinite dimensional twisted local system with a lattice of 
punctures in the d = 4 theory, so that in both cases at least one of the punc- 
tures of such systems lives at the distinguished point of the based loop that 
enters the loop equation. 

A conformal map of the region inside the based loop to a cuspidal fundamen- 
tal domain. 

A regularization of the loop equation at the cusps, that amounts in fact to a 
compactification at infinity of the cusps. 

A reduction of the curvature of the local system at the punctures to an 
Abelian sub-algebra in d — 2 and to a Borel sub-algebra in d — 4, that 
follows by mapping conformally all the cusps (but one) along the same line 
on the boundary of the upper-half plane, by the choice of an axial gauge in 
a direction orthogonal to the boundary of the upper-half plane and by the 
compactification of the cusps at infinity. 

It may be interesting to observe that the construction of the critical equation 
for the quantum holographic effective action, which implies the holographic 
loop equation, involves considering loops with a marked point taken away 
and compactified at infinity after a change of the conformal structure. From 
this point of view, the critical equation is the boundary theory defined on 
the zero dimensional boundary of the loop (i.e. the marked point), while 
the loop equation defines the theory on the bulk. The definition and the 
meaning of all these ingredients will be explained in the following sections. 
In section 2 we recall the holographic form of the loop equation following 
[9]. This section is needed to make this paper self-contained. To check this 
form of the loop equation we reproduce from it to lowest order in pertur- 
bation theory the propagator in the d = 2 case. In section 3 we apply our 
algebro-geometric techniques to the two dimensional case, in order to derive 
from the loop equation for self-avoiding loops the distribution of the eigen- 
values that has been mentioned as the motivation of this paper. In section 



3 



4 we adapt the techniques of section 3 to the four dimensional case and we 
compute the classical and quantum holographic effective action in terms of 
functional determinants. In section 5 we compute f3 , the first coefficient of 
the beta function, from the four dimensional classical holographic effective 
action, finding exact agreement with the perturbative result. In section 6 
we compute the contribution to (3\, the second coefficient of the beta func- 
tion, that arises from the rescaling anomaly which occurs in passing from 
the Wilsonian to the canonically normalised form of the holographic effec- 
tive action. This computation has some interest in itself and it is largely 
independent of the entire holographic construction. It is also perhaps related 
to a question about the beta function of QCD raised in [16]. In section 7 we 
collect some miscellaneous observations, referring especially to analogies in 
the existing literature. In section 8 we state our conclusions. 



We can summarise the basic philosophy in [9] as follows. The loop equation 
in its usual form follows from the observation that the integral of a derivative 
vanishes: 



where the sum over the indices /i, v in the action and over the index a is 
understood. Here A a are Hermitian generators of the Lie algebra and 



ty(x, x;A) is the monodromy matrix of the connection along the closed 
loop C(x, x) based at the point x, i.e. ^(x, x; A) is the Wilson loop. For the 
group SU(N) using the identity: 



2 The holographic loop equation 
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we get: 



/JS f JM 
ZM M exp(-— J TrF^d 4 x)(Tr(-D^ l/ (z)*(x,x;A)) + 

+i f dy v 8W(z-y)(Tr(*(x,y-,A))Tr(V(y,x-,A)) + 

JC(x,x) 

-±Tr(*(x,y;A)*(y,x;A)))) (5) 

where the last term vanishes in the large- ./V limit. The first term is the 
classical contribution to the loop equation, while the second term is the 
quantum contribution. From a modern point of view it is convenient to 
rephrase the loop equation into an algebraic language [17], since in this way 
a more powerful view of what the problem is and of its difficulty is obtained. 
Using the factorisation of gauge invariant operators [18] in the large- iV limit 
and noticing that the expectation value can be combined with the matrix 
trace to define a new generalised trace r, our problem is to find a (unique) 
solution A^(x) to 

= \t{D ii F iu ,(z)*(x,x; A)) + 
+i I dy^ d \z-y)r^(x,y;A))r^(y,x;A)) (6) 

JC(x,x) 

for every closed contour C, with values in a certain operator algebra with 
normalised (r(l) = 1) trace r [19, 7, 20]. Such solution is named the master 
field [21]. Thus the trace r acts on a type II\ von Neumann algebra [17] 
generated by the monodromy operator ^(x, x; A) of Wilson loops based at 
x. This is the loop algebra, that is a representation of the homology algebra of 
based loops at x. Unfortunately, even in the case in which A^x) is Gaussian, 
this algebra is not hyper-finite at iV = oo, that is, it is not the limit of a 
sequence of finite dimensional matrix algebras, being a Cuntz algebra with an 
infinite number of generators [22, 23, 24], which is algebraically isomorphic 
to a free group factor with an infinite number of generators [6]. In the non- 
Gaussian case [19, 7, 20] there is no reason for which things should be easier 
as to the hyper-finiteness property. The point of view pursued in this paper 
is that we should abandon, at first stage, the idea of solving the loop equation 
defined over the entire loop algebra. We should rather solve the following 
problem, whose solution still conveys a lot of physical information, but that 
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is algebraically considerably simpler. For any fixed (self-avoiding, planar) 
based loop, C(x,x), we consider the von Neumann algebra generated by 
ty(x, x; A) only, which we want to determine from the loop equation restricted 
to C(x,x). This is a commutative von Neumann algebra, for which, as it is 
well known, it there exists a structure theorem that is equivalent to measure 
theory plus spectral theory of self-adjoint operators [25]. In particular every 
commutative von Neumann algebra is type / and thus hyper- finite [26] . The 
trace on such algebra is a measure determined by the distribution of the 
eigenvalues pc(^), counting multiplicity [25]: 



Thus our complicated algebraic problem reduces to finding the distribution 
of eigenvalues for any given fixed (self- avoiding, planar) based loop C(x,x). 
It should perhaps be repeated that now both the trace and the (commuta- 
tive) algebra are hyper-finite, so that the trace r, contrary to the original 
problem, is completely known as the large- N limit of the normalised finite 
dimensional matrix trace j^Tr^- In this respect we should mention that, al- 
ready many years ago, the Migdal-Makeenko equation for any given loop was 
written in terms of the distribution of the eigenvalues in a way that makes 
(implicit) use of the algebraic commutative structure associated to iterating 
a given loop [8] . For self-avoiding loops in d = 2 in the weak coupling phase 
on a sphere and for small g 2 A the explicit answer for the distribution of the 
eigenvalues is the formula mentioned in the introduction. In this paper we 
develop methods to solve for the distribution of the eigenvalues in d — 2 
and d = 4 directly from the loop equation for planar self-avoiding loops. As 
mentioned in the introduction our basic idea is holography: we would like 
to map holographically the boundary loop equation into an equivalent bulk 
holographic effective action, whose critical equation determines the distribu- 
tion of eigenvalues, pc, via the eigenvalues of the curvature of a (twisted) 
local system. It may be guessed that in doing so the Cauchy theorem will 
play a key role, being a case (rather spectacular) of planar holography "ante 
litteram". We manage to change variables in the loop equation in such a 
way that in the new variables the quantum contribution be a regularized 
residue, loop independent for self-avoiding loops. This has been achieved 
in [9], changing variables in such a way that functional differentiation of 
the monodromy in the new variables, that is the basic operation to get the 




(7) 
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loop equation, produces the contour integral of a Cauchy kernel instead of a 
delta-like contact term. This makes us a first step closer to find an equivalent 
holographic quantum effective action. 

The second step, by the way, following the d — 2 analogy of a critical equa- 
tion that determines the distribution of eigenvalues, is to make the quantum 
contribution vanishing in the large- N limit. This is achieved in this paper by 
a suitable conformal mapping: the key point is that evaluating the regular- 
ized residue does not commute with the conformal mapping. For our aim, we 
start with representing the partition function as an integral over microcanon- 
ical strata, introducing a suitable resolution of identity [9] . These strata are 
characterised by given levels of the curvature of the gauge connection. Then 
we change variables to a holomorphic gauge in which the curvature F ac- 
quires the form F — d(...). The new holographic loop equation follows. The 
last step, to annihilate the quantum term, is a conformal map to a cuspidal 
fundamental domain. We now describe our procedure in more detail, pro- 
ceeding in parallel in d = 2 and d = 4. The d — 4 case is not substantially 
more complicated than the d = 2 case from a purely holographic point of 
view. As a first step we would like to change variable in the functional in- 
tegral from the gauge connection to the curvature. Formally this is done by 
means of the resolution of identity: 

1 = J Dfi5(F A - n) (8) 

in d — 2 and: 

1 = J D^J(F^(A) - /v) (9) 

in d — 4, that we prefer to write decomposing the curvature into its anti- 
selfdual (ASD) and self-dual (SD) parts: 

i = Jd^j(f- v (a)-^)x 

x / D^6(F+(A)-vl+,) (fO) 

In [9] the resolution of identity associated to the ASD constraint only was 
imposed for reasons that will be cleared in section 4. The SD and ASD 
constraints are written in two dimensional language as Hitchin equations 
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[27, 28]. In particular the ASD constraint is interpreted as an equation for 
the curvature of the non- Hermit ian connection B = A + D = (A z + D u )dz + 
(A z + D a )dz and a harmonic condition for the Higgs field \& = — iD = 
—i(D u dz + D u dz). A is the projection of the four dimensional Hermitian 
connection onto the (z = Xq + ixi, z = x — ix{) plane of the planar loop 
and D is the projection of the four dimensional ant i- Hermit ian covariant 
derivative onto the orthogonal (u = x 2 +ix 3 , u = x 2 —ix 3 ) plane. In this paper 
we choose the following notation as far as the complex basis of differentials 
dz = dx + idx\ and derivatives d = ^ = \(-^ — is concerned. Thus, 
for example, A z = \{Aq — iA\). We should notice that the observable in 
our loop equation is somehow adapted to the microcanonical resolution of 
identity and thus in d — 4 the Wilson loop involving the non-Hermitian 
connection, B, is considered. Yet, since at the end we would like to compute 
planar Wilson loops for the Hermitian connection A, a more general kind of 
observables is needed. Indeed in section 4 we will consider B x = A + XD in 
the limit A — > 0, which we refer to as the unitary limit. In this section for 
simplicity we limit ourselves to the first case. We assume a partial Eguchi- 
Kawai reduction from four [29, 30, 31, 32, 33, 34] to two dimensions [9], 
that implies a rescaling of the classical action by a factor of N^ 1 [34]. This 
factor is needed because the partial Eguchi-Kawai reduction re-absorbs some 
space-time degrees of freedom into the colour degrees of freedom. In our case 
a two dimensional torus is reduced to a point. If the torus is commutative, 
N 2 is given by -^p- / d 2 xd 2 p = In the non-commutative case, instead, 

N 2 = Tr(l) = J2n<N 2 1- The loop equation for B, in the partially reduced 
d — 4 theory, reads: 

= j DBaeM .E.S YM){ Tr(^^ { ^B)) + 

-if dy a S^(z-y)TrCi(x,v,B))TrCi( y ,x;B))) (11) 

JC(x,x) 

with a = 1, 2, so that for these variables it is as difficult to solve as for the 
original four dimensional connection A^. After implementing in the func- 
tional integral the resolution of identity by means of the gauge orbits of the 
microcanonical ensemble mentioned before, we change variable in the loop 
equation from the connection to the corresponding curvature in the holomor- 
phic gauge. Thus we choose for the connection A = a z dz + dzdz in d = 2 and 
B = b z dz + b z dz in d — 4 the gauge d z =0 and b z =0 respectively, performing 
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a gauge transformation in the complexification of the gauge group. This last 
change of variable is made in order to compute explicitly the quantum con- 
tribution as a (regularized) residue. It is well known by the Cauchy formula 
that the line integral along a closed loop of a holomorphic function times the 
Cauchy kernel with pole not lying on the loop depends only on the winding 
number of the loop around the pole and on the value of the holomorphic 
function at the pole of the Cauchy kernel: 

f(z)Ind c (z) = -L / ^-dw (12) 
2m Jc z — w 

where Indc{z) is the winding number of C around z. The change of vari- 
ables to the holomorphic gauge implies that the Cauchy kernel is generated 
by functionally differentiating the monodromy of the connection in the loop 
equation with respect to the integration variables in the functional integral. 
However, even if we reduce the computation of the quantum contribution to 
the evaluation of a line integral of a Cauchy kernel, there are two obstructions 
for the residue theorem to apply. The first one is that the monodromy is not 
a holomorphic function of its endpoints. The second one is that gauge invari- 
ance requires that the functional derivative with respect to the integration 
variable be taken at a point of the loop, for the result to be non-vanishing, i.e 
at a singularity point of the Cauchy kernel. In fact we are going to compute a 
regularized residue obtained as the line integral of a distribution. This solves 
the regularity problem and the holomorphic problem at the same time. In 
d = 4 for B = A + W the microcanonical representation holds in the form: 

r N 
Z= 5{F B - fi)5(<r A V - v) exp(-—S Y M)DBDfj,Dv (13) 
J t±g 

where we have re-casted the ASD microcanonical resolution of identity in 
the Hitchin form of a curvature equation for B plus a harmonic condition for 
the Higgs field for later convenience. The curvature, Fg, is not Hermitian 
in general and thus we may include in the functional integral the resolution of 
identity for its adjoint, in such a way that the integration measure is DfxDp 
instead of D\x. Yet, since the monodromy of B is in fact a functional of \i 
only, omitting Dp does not affect the loop equation for B. Our convention is 
that \i = [i + n — n is the curvature of B in the basis dx Adxi with v = n + n. 
Thus /i° is Hermitian in this basis. Only if we choose a basis different from 



9 



the real one, indices are added to fx to make the choice of basis explicit. From 
now on in this section we mention only the d = 4 case, since the d = 2 case 
is obtained trivially substituting the connection B with the connection A in 
d = 2. We now change variables to the holomorphic gauge, in which the 
curvature of B is given by the field //, obtained from the equation: 

F B - n = (14) 

by means of a complexified gauge transformation G(x; B) that puts the con- 
nection B = b + b in the holomorphic gauge b = 0: 

Bb z = -i^ (15) 

where // = GfxG -1 . The partition function becomes: 

r N Da 

Z = J S(F B - fi)5(d* A ^ - v) eM- — S YM )DB^D f i'Dv (16) 

The integral over B can now be performed and the resulting functional de- 
terminants, together with the Jacobian of the change of variables to the 
holomorphic gauge, absorbed into the definition of T. F plays here the role 
of a "classical" action, since we must integrate still over the field //. We 
may call T the "classical" holographic action, as opposed to the quantum 
holographic effective action, to be found in section 4. T is written explicitly 
in section 4 in terms of functional determinants. The partition function is 
now: 

Z = J exp(-r)D// (17) 

To realize our aim of getting the Cauchy kernel it is convenient to study 
the loop equation for the Wilson loop involving the connection b, thought 
as a functional of B corresponding to gauge transforming B into the gauge 
6 = 0. Such a gauge transformation belongs to the complexification of the 
gauge group and it is rather a change of variable than a proper gauge trans- 
formation. However, because of the property of the trace, for closed loops, it 
preserves the trace of the monodromy. This allows us to transform the loop 
equation thus obtained into an equation for the monodromy of B. In our 
derivation of the loop equation, a crucial role is played by the condition that 
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the expectation value of an open loop vanishes. In [9] two slightly different 
ways of achieving the vanishing of the expectation value of open b loops were 
presented. We may thus derive our loop equation: 



The only non-trivial case is when w lies on the loop C. In this case the loop 
equation can be transformed easily into an equation for B. It is clear that the 
contour integration in the quantum term of the loop equation includes the 
pole of the Cauchy kernel. We need therefore a gauge invariant regularization. 
We proposed in [9] several slightly different ways of regularising the Cauchy 
kernel, that are essentially equivalent concerning the holographic form of 
the loop equation. The first one consists in analytically continuing the loop 
equation from Euclidian to Minkowskian space-time. Thus z — > i{x + + it). 
This regularization has the great virtue of being manifestly gauge invariant. 
In addition this regularization is not loop dependent. A natural, but loop 
dependent, regularization of the quantum contribution to the loop equation 
is a ie regularization of the Cauchy kernel in a direction normal to the loop. 
We have the two possibilities of taking the internal or the external normal. 
We check at the end of this section that the solution of the loop equation 





(18) 



that in the large- N limit reduces to: 




(19) 
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does not depend indeed on this choice to lowest order in perturbation theory 
in the d = 2 theory. Alternatively, as suggested in [9], we may perform a 
conformal mapping of the region encircled by the Wilson loop to the upper- 
half plane, followed by regularization by means of analytic continuation to 
Minkowskian space-time or ie regularization in a direction normal to the 
loop. A more sophisticated form of this regularization will be presented in 
the next section, where we will map conformally the region whose boundary 
is the loop with a marked point to a cuspidal fundamental domain. Quite 
interestingly we will find that the regularization of the loop equation does 
not commute with this conformal mapping. Indeed we will take advantage 
of this fact to annihilate the quantum term, after having accounted for the 
effects of the conformal map everywhere else. In any case the result of the 
ie regularization of the Cauchy kernel is the sum of two distributions, the 
principal part plus a one dimensional delta function: 

-d-\w x -y x + ie) = {2k)~\P{w x - y^ 1 - mS{w x - y x )) (20) 
The loop equation thus regularized looks like: 

= J D^exp(-r)(Tr(^|— *(x,x;b)) + 
- / dy x (2ir)~ 1 (P(w x - y^' 1 - m8(w x - y x )) x 

J C(x,x) 

xTr(^(x,y;b))Tr(^(y,x;b))) (21) 

Being supported on open loops the principal part does not contribute and 
the loop equation reduces to: 

= J D f i'exp(-r)(Tr(j^*(x,x;b)) + 
+ / dyj-5(w x - y x )Tr^(x, y; b))Tr(V(y, x; b))) (22) 

JC(x,x) Z 

Taking w — x and using the transformation properties of the b monodromy 
and of fx(x)', the preceding equation can be rewritten in terms of the con- 
nection, B, and the curvature, /i: 

r 5T 
= J D^'exp(-T)(Tr(j-^^(x,x-,B)) + 

+ / dy x % -5{x x - y x )Tr(V(x, y; B))Tr(*(y, x; B))) (23) 

JC(x,x) Z 
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where we have used the condition that the trace of open loops vanishes to 
substitute the b monodromy with the B monodromy. This is our final form 
of the regularized Euclidian loop equation (there is an analogous form in 
Minkowskian space-time). Let us notice the sign ambiguity in the quantum 
contribution that depends on the choice of the sign of it. Because of the 
product of traces, in general the quantum contribution does not have the 
same operator structure as the classical term. In addition the quantum con- 
tribution is loop dependent in general. However, for self-avoiding loops, the 
quantum contribution is just a linear topological term added to the "classi- 
cal" action, T, provided the loop equation is interpreted in a strong sense, 
that is as: 



An interpretation of this kind was attempted in [9], where the central term, 
being topological, was thought to be essentially irrelevant or cancelled against 
an anomalous phase in the effective action. Unfortunately, following the in- 
terpretation of [9], we were unable to reproduce to lowest order in pertur- 
bation theory the gauge propagator. The reason is that the central term 
is absolutely essential to obtain the gauge propagator to lowest order from 
the loop equation for self-avoiding loops, if the loop equation is interpreted 
weakly as: 



In fact the role of the central term in getting the correct propagator is en- 
tirely analogous to the role of the contact term in the original loop equation 
(Eq.(5)). Hence we might conclude that we have come to a loose end despite 
our sophisticated changes of variable: perhaps we have simply rewritten our 
loop equation in an exotic way for exotic variables. In the two next sections 
we will see that, thanks to some new ingredient, this is not in fact the case. 
We end this section checking to lowest order that the gauge propagator does 
follow from the loop equation in the holographic form in d = 2, when it 
is interpreted in a weak sense. For computational convenience we con- 
sider the holographic loop equation in d — 2 followed by Eguchi-Kawai 
reduction. We interpret it in a weak sense and we expand in powers of 
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(24) 




(25) 
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fj,(x) = exp(ipx)ji(0) exp(-ipx). We get for a small loop: 




where the factor of A r 2 ~ 1 = N' 1 takes into account the Eguchi-Kawai reduc- 
tion from d — 2 to d — and L 2 is the area of the space-time torus. A 
somehow intriguing feature of our computation, based on the Eguchi-Kawai 
reduction, is that we perform it in a way that implies (to lowest order) a 
hyper-finite trace. Yet this does not imply necessarily hyper-finiteness of 
the algebra, that would follow instead from uniform hyper-finiteness of the 
trace [26]. Indeed even for the solution in terms of free random variables in 
d = 2, that certainly generate a non-hyperfinite algebra, a solution admitting 
a hyper-finite trace has been proposed [7]. Thus if we set: 



for i > k with fj,^ = and /j, = fi*. Since from Eq.(8) it follows, to lowest 
order in powers of fi, in the gauge D a 5A a = 0: 



r = lim — Trjv 



(27) 



we get for fi from Eq.(26) the normalisation: 



tok(x) = -exp(i(pi -pk)x) 



(28) 



A z = id(-A) 
A- z = -id(-A) 



i 



i 



(29) 



we get for the propagator to lowest order: 



r(A z (x)A- z (y)) 




(Pi ~ Pk) 



1 



-exp(i(pi-p k )(x-y))) 




(30) 



as it should be in the quenched theory. 
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3 The two dimensional case: reduction to an 
Abelian sub-algebra 



The aim of this section is to solve the holographic loop equation of the previ- 
ous section by means of a certain Abelianization map in the two dimensional 
case. The advantage of considering the two dimensional case first is that in 
two dimensions we already know the exact answer, so that we can test our 
ideas about solving the loop equation for self-avoiding loops. As mentioned 
in the introduction, to realize the Abelianization map we need: 
a local system; 

a uniformization to a cuspidal fundamental domain; 

a regularization at the boundary cusps that amounts to a compactification 

of the cusps at infinity; 

a clever choice of the gauge. 

The local system was introduced in [9] basically to have a mathematically 
well defined model of the moduli space of the master field. Yet, it turns out 
that the need of a local system and in particular of the associated singular- 
ities has a deeper meaning. The local system furnishes a lattice or adelic 
interpretation of the microcanonical localisation: 



From the point of view of the functional integration the curvature of the 
connections associated to local systems is dense in the sense of distributions 
in the space of curvatures. For local systems the loop equation acquires the 
following form: 



F A = ^2^ P S 2 (X - X p ) 



(31) 



V 



in the holographic loop equation: 




(32) 



= J J| Dfi' q exp(— r)(Tr(- — ^(x p , x p ; A)) + 




dy z -d-\xp - y)Tr(V(x p , y; A))Tr(*(y, x p , A))) (33) 
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where now it is necessary that the marked point of the based loop coincides 
with a puncture to get a non-trivial loop equation. It can be checked by direct 
computation that integrating in the functional integral over non-Abelian local 
systems leads the same result for the Wilson loop as in the lattice theory: 

t(V(A;x,x)) = Z- 1 f HdgpTrigpeMie^g^ 1 ) x 
J P 

x «p(- E ^L T < e l)) II II I «p(4) - ™?H)\de P (34) 

p Z 9 a p ijtj 

where the product over p is restricted to the lattice points internal to the 
loop and a~ 2 = ^2 = 5^(0). The dg p integrals can be easily performed 
leading to the correct lattice result [3] (the only difference is the Wilsonian 
lattice action instead of its formal continuum limit) : 

r N 
t(*(A;x,x)) = Z l \{ /(Tr(exp(^ p ))exp(-- TT Tr(^)) x 
p J zg a 

x Yl I exp(i9p) - exp(i8 3 p )\d9 p 

that is equivalent for r(\l/) to the formula for the distribution of eigenvalues 
in the introduction for small g 2 A and in the scaling limit Nq — > 00 with 
Ncci 2 = A = constant, where Nq is the number of punctures inside the loop 
C. We will see now how this result is reproduced by the Abelianization map of 
the holographic loop equation. The key point is that identifying the marked 
point of the based loop with a puncture contains implicitly the possibility of 
a change of the conformal structure. Formally the Wilson loop is invariant 
under re-parameterisation of the boundary that can be extended to conformal 
transformations of the region encircled by the loop, but in fact we will change 
the conformal structure in a way that is equivalent to attach to the loop, in 
a neighbourhood of x, an infinitesimal strip going to infinity. This does 
not alter the Wilson loop because of the zig-zag symmetry [35]. Of course 
the classical action in d = 2 is not conformally invariant, so that if we are 
going to compute the classical action in terms of a conformally transformed 
local system we must transform the classical action properly. Coming back 
to the loop equation for local systems, it is of the utmost importance that 
the regularized residue with the new conformal structure vanishes in the 
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loop equation. We now explain why, first proceeding heuristically and then 
by direct computation. Our lattice gives rise to a punctured sphere with 
a based contour, C(x,x), that determines an internal region (the smaller 
one on a large sphere), fl x , with the topology of a disk with some punctures 
inside and at least one on the boundary, that coincides with the distinguished 
base point of the loop, x. Q x can be mapped conformally to a cuspidal 
fundamental domain on the upper-half plane with all the cusps but one, for 
example the cusp on the boundary, on the y = axis and the remaining 
one at y = oo. On the other hand also the region external to the loop can 
be mapped conformally in a similar way. Heuristically this conformal map 
allows us to diagonalise the curvature at the cusps, thanks to a residual gauge 
symmetry (see below), provided the cusps are compactified. In general the 
compactification is not possible, otherwise every local system on a sphere 
would be Abelian. However, in our case, the compactification is required 
since the marked point, x, in origin belongs to the closed loop C(x,x) and 
thus it is not a puncture. In addition, assuming rotational invariance on 
the sphere, if one cusp is compactified also all the remaining ones should 
be. On the fundamental domain all the cusps but one are on the x-axis and 
hence, choosing an axial gauge along the y-axis, there is a residual gauge 
symmetry of making gauge transformations along the to diagonalise 

the curvature at the cusps on the In addition, since our local systems 

define a representation of the fundamental group of a punctured sphere, also 
the curvature at the remaining cusp at y = oo must be diagonal. Thus we 
get an Abelian local system. But Abelianization can only be consistent with 
a vanishing of the regularized residue in the loop equation, since an Abelian 
system becomes classical in the large- N limit. Thus there is no quantum 
term in the large- N loop equation on the cuspidal fundamental domain, and 
being gauge invariant, the loop equation is to coincide before and after gauge 
fixing. Hence there never is a quantum residue in any gauge. Notice that the 
Abelianization could not be performed before the conformal mapping, since it 
is essential that the cusps be on the same line and compactified. Now we must 
check directly that the quantum residue vanishes in any gauge. The cuspidal 
domain is a polygon defined by the uniformization theory of Riemann surfaces 
with punctures and boundaries. A unified approach in which the boundaries 
are treated in a way similar to the punctures has been given by Penner [36]. 
Remarkably Penner approach involves the choice of at least one puncture on 
the boundaries, precisely as required by our interpretation of the holographic 
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loop equation in terms of local systems. The boundary arcs of the polygon 
that uniformizes our region Q x are oriented in the following way. The couple 
of arcs that end into a cusp corresponding to the punctures in the interior 
of Q x have opposite orientations on the fundamental domain, since they are 
in fact identified by the gluing map that reconstructs from the polygon the 
punctured Riemann surface: this identification creates tubes out of such 
cusps. As a consequence these couple of arcs share the same orientation on 
the reconstructed Riemann surface. On the contrary, the arcs ending into 
the cusp lying on the loop share the same orientation on the polygon since 
they are not glued together: these arcs are associated to an infinitesimal 
strip going to infinity. As a consequence these couple of arcs have opposite 
orientation on the strip to infinity. This difference in orientation plays a 
crucial role in evaluating the regularized residue at the internal cusps and 
at the boundary cusps. In fact in the first case the regularized residues 
associated to the two asymptotes of the cusp sum up to 1 because of the 
same orientation of the asymptotes on the Riemann surface. In the second 
case the sum is 0, because the opposite orientation of the asymptotes on the 
Riemann surface. We may look at the last fact as just another consequence 
of the zig-zag symmetry. Thus the quantum contribution vanishes. We can 
summarise our argument as follows. The local system is invariant under the 
conformal map, while the classical action transforms in a definite way. The 
marked point of the based loop becomes a cusp of the fundamental domain. 
This cusp is obtained adding an infinitesimal strip going to infinity to the loop 
in a neighbourhood of the marked point. Because of the zig-zag symmetry, 
at this cusp the quantum regularized residue vanishes. We can take into 
account the preceding argument to derive a new form of the holographic 
loop equation that leads to a critical equation for a holographic quantum 
effective action. The resolution of the identity involves now local systems on 
the cuspidal fundamental domain with coordinates (t,t): 



where the superscript ^ refers to fields defined on the fundamental domain. 
It is this new resolution that is inserted into the functional integral: 




(35) 



= 



TrF 2 A d 2 x)5(Ff - ^)D^ 



(36) 
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Then all the steps go through as in the second section. The universal cover 
of the fundamental domain is the upper-half plane U. On U we choose the 
gauge A y = 0, that leaves, as a residual gauge symmetry, gauge transfor- 
mations that are y independent. In the gauge A y = the determinant due 
to localisation, i.e. the one obtained integrating with respect to the gauge 
connection, A, the delta functional in Eq.(36), and the Faddeev-Popov de- 
terminant are both field independent and cancel each other. We can use 
the residual gauge symmetry to fix the gauge /ijf = (the label ch means 
the non-diagonal part) at the cusps on the x-axis. The associated extra 
Faddeev-Popov determinant is the square of the Vandermonde determinant 
of the eigenvalues of the curvature of the local system at the punctures. Then 
the holographic quantum effective action reduces to: 



where the ^ superscript refers to the domain of definition of the lattice field, 
hp, and we have set a = with a the lattice spacing corresponding to the 
cutoff A of the theory, that comes from the product of delta functions at 
the same point in the classical action. The last term is the logarithm of the 
Jacobian to the holomorphic gauge of the Abelian local system. However it 
vanishes identically in a gauge in which the curvature is Abelian (see at the 
end of next section). It should be noticed that T q is expressed as a functional 
of the local system on the fundamental domain. This involves in the classical 
term a change of the metric since the classical action is not conformally 
invariant. Since on the fundamental domain the quantum term vanishes, the 
loop equation for self-avoiding loops reduces to: 



N 



E E l£(p)W )a -EE WW - hf>\ + 

-logDeti^) 




2g 2 a 2 



(37) 




^{x p ,x p ,A)) 



(38) 



Eq.(38) is implied by the critical equation: 



o 



(39) 




19 



Getting a saddle-point morally implies that we have somehow integrated 
away the order of iV 2 non-Abelian degrees of freedom to obtain the effective 
action for the order of iV remaining eigenvalues. It is the conformal map that 
allowed us to diagonalise the curvature at the cusps. But this works only for 
a system in which the curvature has delta-like singularities. Hence local sys- 
tems know about quantum field theory ! It is perhaps this the very reason for 
the occurrence of Hitchin systems in d — 4 quantum field theories (see section 
4 and section 7). Notice that the Abelianization works only for self-avoiding 
loops. In case of self-intersection the marked point would be ramified and 
globally there would not be enough residual symmetry in an axial gauge on a 
ramified covering to diagonalise all the cusps. On the other hand the regular- 
ized residue in the loop equation does not vanish in general at a ramification 
point. It remains to be seen that our critical equation coincides with the 
large- N saddle point in Eq.(l). To do so we need a better understanding 
of the uniformization map. It turns out that the relevant mathematics is 
still the one of the moduli space of bordered Riemann surfaces following [36]. 
As already noticed, Penner version of the Teichmuller theory of bordered 
surfaces involves the choice of at least a puncture and a horocycle arc in a 
neighbourhood of the puncture on each boundary of the surface, in analogy 
with the theory of punctured surfaces that involves a horocycle around each 
internal puncture. More generally we may consider the case in which the 
loop in the loop equation intersects a number of punctures of our lattice. It 
is clear that the preceding arguments about the vanishing of the quantum 
residue apply in this more general situation, since they depend indeed on 
the local structure around each puncture. But now, since all the punctures 
on the loop are attached to infinitesimal strips going to infinity, the circle 
at infinity contains a lattice (or adelic) image of the original loop through 
the conformal map. Quadratic differentials [37] can be used to construct the 
uniformization map from the generic region that occurs in the loop equation 
to the cuspidal fundamental domain. The basic relation between quadratic 
differentials, q, and the uniformization map is: 

| - ^ (*) 
We need therefore the standard form of a quadratic differential near a cusp: 

* 1 (41) 



dz 2%iz 
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where L is the length of the horocycle arc around the cusp. Since this expres- 
sion is infinite at the cusps it must be regularized and suitably interpreted. 
In particular it depends crucially on what the cutoff is on the fundamental 
domain near the cusps. This may be difficult to understand in general, but it 
seems easier in the case of a circular loop. In this case we have essentially a 
circle that is mapped (adelically) into the circle at infinity. This is a cylinder, 
i.e. a punctured disk, which is mapped by the uniformization map to a strip 
in the upper-half plane. In this case the uniformization map is: 

1 = ^n° 9{z) (42) 

Thus we get: 

\^m 2 = ^ = — 2 - N c (43) 

oz a 2 rra 2 

where R is the radius and A the area of the disk, while a is the radius of 
a little disk around the puncture. Thus Nc is the number of lattice points 
inside the disk. When Eq.(43) is inserted into Eq.(37) the correct distribution 
of eigenvalues, given by the saddle point equation for the effective action in 
Eq.(l), is obtained after noticing that 9 l in Eq.(l) is related to h l p = h l (by 
rotational invariance on the sphere) in Eq.(37) by: 

Q i = Nch 1 (44) 

since the magnetic flux through the loop counts the number, Nc, of punctures 
inside. In fact the holographic action in terms of 6 l becomes Z)p( 2g 2jv ca 2 ^* & l2 ~ 

log\6 l — Qi\ + constant) which differs by an irrelevant overall factor and 
an additive constant from the one in Eq.(l). In this paper we have ignored 
the contribution of the region external to the loop and the corresponding 
conformal map. We will consider it elsewhere. It should be noticed that the 
eigenvalues of the curvature at the punctures of the local system are in fact 
defined modulo 2-n. We can take into account the periodicity of h p summing 
over appropriate windings. This is expected to lead to the known phase tran- 
sition for g 2 A large in the large iV limit [3, 4, 5], but in that case the external 
region must be considered too. 
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4 The four dimensional case: reduction to a 
Borel sub-algebra 



The key points in the previous section can be extended to the planar four 
dimensional case. We list here the needed modifications. We should re- 
mind that in our approach the observable in the loop equation is adapted 
to the microcanonical resolution of identity. In the four dimensional case 
the resolution of identity involves the SD or ASD curvature. Thus, for ex- 
ample in the ASD sector, the connection that enters the loop equation is 
B = A + D, which is non-unitary. However for physical reasons we are in- 
terested to compute the Wilson loop for the unitary connection A. Thus we 
introduce B x = A + (\D u dz + \D u dz) where A is a section (possibly con- 
stant) of a holomorphic line bundle and we take the limit A — ■> 0. We refer 
to this limit as the limit of unitary Wilson loop or the unitary limit in short. 
In the four dimensional case it is most convenient to take the space-time to 
be a product of a two dimensional sphere by a non-commutative torus in the 
limit of infinite non-commutativity [9]. This limit is known to be equivalent 
to the usual commutative theory in the large- N limit [31]. Although several 
different space-times may be considered, in this paper we make this choice 
because the associated SD or ASD equations look like a kind of infinite di- 
mensional vortex equations (see below). Heuristically we find appealing the 
occurrence of vortices for the hope of reproducing both the area law at large 
distances and the Coulomb law at short distances. The choice of the sphere 
is also modelled on the analogy with the two dimensional case, especially for 
simplifications due to the non-existence of non-trivial cycles, but for the ones 
associated to the punctures. Indeed on a sphere all the moduli of the local 
system are the local moduli. When we take the unitary limit in our loop 
equation, some care is necessary, since the limit has to be taken in such a 
way that the correct four dimensional information survives, for example in 
the beta function. The first coefficient of the beta function (see next section) 
can be computed in two different but related ways, that we now explain. 
Since the theory lives on a product of a sphere by a non-commutative torus, 
the curvature equation involves a central term, H, equal to the inverse of 
the parameter of non-commutativity, 9. This occurs because, once the gauge 
connection is required to vanish at infinity up to gauge equivalence, the only 
term that survives in the curvature at infinity is the commutator of the 
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derivatives on the non-commutative torus, that is H . In turn H vanishes as 
j; in the large-iV limit. Therefore, in the case of B x , our centrally extended 
and A twisted ASD curvature equation reads: 

F A - t\X\ 2 ^ 2 = £ n°/ 2 \x - x p ) + HI 
v 

xBa^p = ^2 ^ n p^ 2 \ x — x v) 

V 

\d A ^j = ]T \n p 5 {2) (x - x p ) (45) 
v 

where we have set D = i^fr, we have rescaled n by a factor of A to ensure 
the finiteness of ip in the A — > limit and analogously in the SD case. The 
central extension H is referred to in this paper as the twist of the local sys- 
tem. The A rescaling instead is a twist of the Hitchin system considered as 
a twisted Higgs bundle. In our interpretation of the preceding equations as 
vortex equations, the central extension H is related to the non-vanishing of 
the Higgs field at infinity, while the zeroes of the Higgs field may arise from 
twisting by the factor of A. In fact this is precisely what we require in the uni- 
tary limit. We choose a A that has a lot of zeroes in a compact set containing 
the loop C, in order to make the Higgs field vanishing small in a neighbour- 
hood of the loop to ensure unitarity of the monodromy along C, and that 
converges to 1 at infinity, to keep the information about the four dimensional 
nature of the theory. As to the beta function, the simplest case is A = 1. In 
this case, if we are interested only in the ultraviolet logarithmic divergences 
for the collective field /x° and not in the critical equation for the holographic 
quantum effective action, we need not to decompose the collective field /i° 
into a sum of delta distributions, that is equivalent to introduce the local 
system, and in fact (3o can be found directly from the ultraviolet divergences 
of the "classical" holographic action, T, looking at the Tr(/j? 2 ) counter-term 
(see next section), solving in perturbation theory, at first order in power of 
the local curvature, the preceding equation around the jjP = n = solution. 
However our observable in the case A = 1 is not physically interesting. In 
addition it is trivial at least at first order in perturbation theory and in the 
large- limit. This occurs because the contribution of the propagator of the 
field A z is exactly cancelled by the one of A u because of the different factors 
of i in the Wilson loop. In fact at A = 1 the Wilson loop is probably topolog- 
ical in the large- N limit. The limit A — > is the most physically interesting. 
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In this limit H must be left untouched at infinity since it is essential to keep 
the correct four dimensional information while n is expected to be rescaled 
by a factor of A. A little thought shows that this may happen if A has a 
zero at each puncture in a compact set containing the loop and converges 
to 1 at infinity as already anticipated. This implies that B x converges to 
(A z + d u )dz+ {Az + d n )dz at infinity and to A on a compact set on the sphere 
in the complement of infinity. In this case it is more convenient the lattice 
or adelic interpretation of the curvature equation in order to obtain (3 Q . It is 
natural to interpret the adelic theory as defining a scaling limit as in usual 
lattice gauge theories. In such theories there are no ultraviolet divergences 
because of the finite lattice spacing, but all the divergences appear as loga- 
rithmic infrared divergences at a scale much larger than the lattice spacing 
but still smaller than the inverse of Aq C d, that goes to infinity in the scaling 
limit, i.e. when the coupling constant goes to zero. If we are interested in the 
infrared logarithmic divergence we should employ in our computations the 
asymptotic value of A at infinity, i.e. 1, and hence our computation reduces 
to the one in the case A = 1 everywhere. For our twisted local system the 
region of asymptotic freedom is thus the large distance region on the sphere. 
This perhaps resembles the holographic UV-IR duality already encountered 
in [10]. From a technical point of view, as far as /3q is concerned, the com- 
putation for A = 1 and the one for A — * are identical once it is observed 
that the in the second case we are actually looking at a neighbourhood of 
infinity. Thus in this paper we report only the first case in the next section. 
We have seen that having introduced non-commutativity leads inevitably to 
a central extension in the curvature equation and thus for A = 1 we get a 
central extension of an infinite dimensional representation of the fundamental 
group of the punctured Riemann surface. This central extension cannot be 
decomposed in general into a U(l) Abelian system and a flat connection as 
in finite dimensions. Thus we may wonder which are the finite dimensional 
approximations of this system. The answer is quiver bundles with parabolic 
singularities [38]. In d — 4 the additional difficulty arises that is not possi- 
ble to impose the two SD and ASD constraints independently at the same 
time for an irreducible connection in the functional integral. Thus either we 
impose only the ASD part of the resolution of identity as we did in [9] (we 
are not finding a universal master field, since our observable is adapted to 
the resolution of identity) or we introduce orbifold models as follows. There 
are orbifold models whose large- N limit is equivalent to the large- iV limit 
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of QCD. We explain what a orbifold is in this context following [39]:" A 
certain parent gauge theory is chosen", in our case pure SU(N) gauge the- 
ory; "the orbifold theory is simply given throwing away all fields that are 
not invariant under a discrete subgroup of the gauge symmetry " in our 
case Z 2 . "The resulting theory has the remarkable property that at large- TV 
its perturbation series is the same as the parent theory, up to some simple 
rescaling of N" . Sometimes, using the loop equation, it is even possible to 
show non-perturbative equivalence between the parent and orbifold theory. 
This is the case for our Z 2 orbifold, that is simply a gauge theory with gauge 
group SU(N) x SU(N). It is clear that if we compute the Wilson loop in the 
(N, N) representation in the Z 2 orbifold theory in the large- N limit, this will 
be the same as the Wilson loop in the adjoint representation of the parent 
SU (N) theory. As we said the necessity of a Z 2 orbifold occurs in d — 4, if we 
want to impose in the functional integral the two SD and ASD constraints 
independently at the same time. We now write the loop equation in terms 
of the lattice field of curvatures in d = 4: 



where a A dependent effective action appropriate for the study of the mon- 
odromy of the operator B x = A + + A0) = A x + i^/ x , with /i x z = 
fi° z2 + \n Z z — Xn Z z, must be introduced. As in d = 2, to get rid of the quan- 
tum term, we must map the lattice loop equation to a cuspidal fundamental 
domain. In doing so we make a conformal rescaling of the Hitchin system as 
dictated by the conformal properties induced requiring the conformal invari- 
ance of the non-unitary monodromy that enters the loop equation: 



C J- 

= J Y[Dfx' q exp(-r)(Tr(—^(x p ,x p ;B)) + 




dy z -d-\x p - y)Tr(V(x p , y; B))Tr(*(y, x p , £?))) (46) 



and analogously in the unitary limit: 





•t. 



(48) 
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Thus, in mapping to the fundamental domain, we are actually rescaling in 
a certain way also H, the central extension, and the orthogonal coordinates 
as well. Now, to get a zero quantum contribution, we must insert in the 
functional integral the Wilson loop and the resolution of identity, expressed 
in terms of the Hitchin system on the fundamental domain (we do not display 
explicitly the dependence on A): 

Z = J 8(Fg> - n (t) - H^l)5((d A ^) w - u {t) ) x 

xeM-^S YM )DB^D^D^ = 



J exp(-r)lV (t) (49) 



As in d = 2 the Wilson loop is invariant by construction, up to the added 
infinitesimal strip to infinity, that is attached and got rid by the zig-zag 
symmetry. The classical action must be expressed in terms of the Hitchin 
system on the fundamental domain as in the d = 2 theory: 



The contribution to Sym m the ASD sector maintains the same form because 
of its invariance under conformal rescaling of the metric. For the SD part 
there might be additional terms due to the fact that a holomorphic twist of 
the Higgs field in the ASD sector is anti-holomorphic in the SD one. The 
holographic quantum effective action reads (we do not display explicitly the 
dependence on A): 

T q = T\ Ay=0 \ FD -J2J2 l °9 Det ( ad ^p)\^-=o ( 51 ) 

The first term, T\ Ay=0 \ FD = ((^S Y M + ^logDet'( y -A A 5 flu + D^D u + iad F - i/ )- 

logDet'jyjj)\A y =o)\FD, is the "classical" holographic action associated to the 
reduced non-commutative theory on the fundamental domain in the axial 
gauge. From a computational point of view it is more convenient to calculate 
the functional determinants in terms of the original Hitchin system before the 
conformal map to the fundamental domain. The two expressions for T should 
coincide up to the conformal anomaly. Thus T| j 4 !/= o|fd = ^\A y =o\sphere + 
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Con formal Anomaly. The structure of =o\s P here will be elucidated below. 
The second term is due to the additional gauge fixing to a Borel sub-algebra, 
u~ = 0, and it is the d = 4 analogue of the Vandermonde determinant of 
the eigenvalues, to which it reduces in the unitary limit A — > 0. The label 
for u p means here lower triangular part, excluding the diagonal, while the 
label + for u p means here upper triangular part including the diagonal. The 
conformal anomaly is the contribution due to the change of metric implicit 
in mapping conformally the region encircled by the punctured Wilson loop 
to the cuspidal fundamental domain. It can be obtained from the exact beta 
function. The holographic loop equation on the fundamental domain reads 
(we do not display explicitly the dependence on A): 

= T(-^fc*(x p ,x p ;B)) (52) 
that is implied by the critical equation: 

5^° (53) 

which we refer to as the master equation. We should notice that, in analogy 
with the two dimensional case, corrections due to the external region and 
to the periodicity of the eigenvalues should be included. We write now the 
"classical" holographic action, T, in terms of functional determinants. For 
the purpose of computing the first coefficient of the beta function it is con- 
siderably more convenient to calculate T in a Lorentz gauge rather than Y q in 
the axial gauge, that is needed to realize the reduction to a Borel sub-algebra 
and, in the unitary limit, to a diagonal one. We expect that the divergences 
of T and T q coincide, since T q is obtained from T just by a conformal map 
and a suitable gauge fixing. Since however the conformal map is singular, 
subtleties might be involved and thus we would like to have a direct under- 
standing of T q . Yet, we will not address this problem in this paper. Part of 
the computation of T in a Lorentz gauge appeared already in [9] . However we 
have computed here more explicitly than in [9] the Jacobian of the change of 
variables to the holomorphic gauge. This Jacobian turns out to be essential 
to reproduce the correct value of fto i n the limit A — > 0, which is the one 
relevant to the physical case of Wilson loops for the unitary connection A. 
In this respect the computation of the beta function that we present here is 
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more relevant to the physics than the one in [9], where the Jacobian to the 
holomorphic gauge was not computed and SD of the master field implicitly 
assumed in the sector of ASD type. In this case it was found in [9] that the 
one-loop perturbative beta function was accounted completely by the local- 
isation determinant (see next section). In the Feynman gauge T is the sum 
of the classical (reduced, non-commutative) Yang-Mills action plus the log- 
arithm of a number of determinants whose origin is as follows. There is the 
determinant that arises from the localisation integral on the micro canonical 
ensemble. The localisation determinant is defined formally as: 

J DA,5{F; U - ) = Det'-\p-d A A) 

= Deif-^((p-d A A)*(p-d A A)) 
= Det'^(-A A 5^ + D„D V + iad F -) (54) 

where P~ is the projector onto the anti-selfdual part of the curvature. The 
' suffix requires projecting away from the determinants the zero modes due 
to gauge invariance, since gauge fixing is not implied in the loop equation, 
though it may be understood if we like to. 

The careful reader may have noticed the unusual spin term, iad F - v , as op- 
posed to liadp^ which arises in the perturbative effective action (see next 
section). The occurrence of F~ v is due to the projector P~ in Eq.(54). The 
coefficient of F~ u is one half of the one of F^ since the delta functional in 
Eq.(54) is approximated by a Gaussian whose quadratic form is the second 
order expansion of (F~ u — //~„) 2 around /i~ u , while the usual background field 
method involves expanding the quadratic form associated to F^ u around the 
background /j,^ (see next section). After inserting the gauge fixing condition 
and the corresponding Faddeev-Popov determinant the localisation determi- 
nant can be defined non-formally in the Feynman gauge as: 

Det'-^(-A A 5^ + D„D V + iad F -) = 
= Urn J DcDA^expi-— J rf 4 xTr(c 2 )) x 
x exp(-l £ J d 4 xTr((F^ - ^ v f )))5{D,5A, - c)A FP (55) 
The result is thus: 

Det'-^i-A^^ + D ii D v + iad F - u ) = Der^{-A A 8„ V + iad F -JA FP (56) 
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where is the anti-selfdual part of the field strength, given by: 

F/xv = — F^ v (57) 

with: 

F*!xv — l^nvapFaP (58) 

In addition the determinant that arises as the Jacobian of the change of 
variables to the holomorphic gauge contributes to the " classical" holographic 
action: 

-logDet'^- f (59) 

which can be computed more explicitly as: 

logDet{\ + [G- 1 — (60) 

where the following equation has been used: 

G^BG = -ib (61) 
The "classical" holographic action, T, reads: 

r = ^S YM +]^ogDet(-^ A 5^ + iad F -J - logA FP - logDet^^ (62) 

where Sym is the Yang-Mills action of the reduced non-commutative theory 
and we assume a corresponding reduction in the other terms where neces- 
sary. The normalisation of T given by the reduced theory is needed to be 
compatible with the two dimensional normalisation of the Vandermonde-like 
determinant in Eq.(51). Notice that the last term in Eq.(62) vanishes if jj, is 
reduced by gauge fixing to an Abelian or to a Borel sub-algebra, since the 
commutator in Eq.(60) involves also the colour indices and the commutator 
of elements in the Borel sub- algebra is nilpotent. 
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5 Beta function: the first coefficient 



We compute in this section the first coefficient of the beta function as it 
follows from the "classical" holographic effective action T. In this section for 
simplicity we evaluate T in the un-reduced theory, rather than in the reduced 
one. The reduction is irrelevant as far as (3$ is concerned, since its compu- 
tation does not involve the second term in Eq.(51). A partial computation 
involving only the part of V composed by the classical Yang-Mills action and 
the localisation determinant was already performed in [9]. Following [9], it 
is convenient to perform the computation in an indirect way, by means of a 
term by term comparison with the usual one-loop perturbative contribution 
to the effective action. For this purpose let us recall the structure of one-loop 
perturbative corrections to the classical action in the Feynman gauge: 

r N r N 

J DcDA, exp(-— J d A xTr{c 2 )) exp(-—S YM )5(D u 5A u - c)A FP = 

N i 

= exp(-—S Y M)Det-z(-A A 5^ + i2ad F jDet(-A A ) (63) 

where we have inserted in the functional integral the gauge-fixing condition 
and the corresponding Faddeev-Popov determinant and, by an abuse of nota- 
tion, we have denoted with A the classical background field in the right hand 
side of Eq.(63). It follows that the perturbative one- loop effective action in 
the Feynman gauge is given by: 

N 1 

^one-loop = ^ S ym + -logDet(- A A S av + i2ad F ^) - logDet{-A A ) (64) 

The perturbative computation of the one-loop beta function [40, 41] is the 
result of two contributions that are independent within logarithmic accuracy 
[42]. The orbital contribution gives rise to diamagnetism and to a positive 
term in the beta function: 

-log(Det-^(-A A 8^)Det(-A A )) = logDet(-A A ) = 

= T2<k* l ° 9( ^J dlxTr(F - f (65) 

where it should be noticed the cancellation of two of the four polarisations 
between the first factor and the Faddev-Popov determinant. The spin con- 
tribution gives rise to paramagnetism and to an overwhelming negative term 
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in the beta function [42]: 
1 



£ Tr{i2ad F ^-/\ A )- l i2ad F ^-^ A y l ) = 



2Tr(iad(fi°)(-A A )-Had(^)(-A A )- 1 ) = 

12 N Mi-) 1 -^ j d'xTriF^f (66) 



12(2tt) 2 



where for later convenience we have expressed the spin contribution in the 
A — > limit in terms of the field /x . Hence the complete result for the 
divergent part of T one _[ oop is: 

from which it follows that: 

From Eq.(62) it can be easily read that the orbital contribution in T one ^i oop 
and T coincide. On the contrary, the spin contribution in T involves only 
the anti-selfdual part of the curvature instead of the curvature itself. Hence 
the spin contribution from the localisation determinant in the unitary limit 
is only one half of the spin contribution in perturbation theory. Thus the 
orbital and spin contributions of the localisation determinant sum up to: 

In fact, remarkably, the missing other one half in the coefficient of the beta 
function is furnished by the Jacobian to the holomorphic gauge in the unitary 
limit, which to second order in yU° contributes to Y: 

logDetil-ifi- 1 — l =fl o,v }) (70) 

where the commutator involves space-time indices and colour as well. From 
the ASD localisation it follows to lowest order in /i: 

^ = _i<9(_A) -1 + <9 U (-A) -1 (71) 
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The second term in the preceding equation does not contribute to second 
order because of rotational and parity invariance in the (u, u) plane. Eq.(70) 
contains traces of powers of a commutator, hence a formal evaluation would 
hardly lead to any logarithmic divergence. For example the contribution 
which we are interested in is the trace of a commutator squared: 

~Tr([(-A)-\ad(^)} 2 ) (72) 

where now the commutator involves only space-time indices. Yet, one term 
is zero in dimensional regularization, being a tadpole: 

Tr{{-A)-\ad{^f) (73) 

while the other term: 

-Tr((-A)- 1 arf( / u°)(-A)- 1 ad( y[ x )) (74) 

furnishes the needed contribution with the correct sign and coefficient to 
combine with the result from localisation exactly into the perturbative (3q. 
This seems to be a non-trivial check of our chain of changes of variable. 



6 Beta function: the contribution of the rescal- 
ing anomaly to the second coefficient 

It is an interesting question as to whether Y reproduces also the second 
coefficient of the beta function. We do not have a definite answer at the 
moment. A contribution to (3i, if it exists, can only come from higher order 
operators, starting with Tr(/x 4 ), that would carry the correct power of g. It 
is then clear that, to relate Tr{ji A ) to Tr(/x 2 ), Eq.(25) should be employed. A 
shortcut may perhaps consist in choosing a // that reproduces to lowest order 
the gauge propagator as in the two dimensional case. A most natural thing 
would be to compare our holographic effective action with the Wilsonian 
effective action of QCD, based on the exact renormalization group [43, 44], 
if we knew it. In principle an ansatz for the exact Wilsonian effective action 
should include all higher order operators as it turns out in our holographic 
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case. Yet the simplest ansatz that includes only the lowest order operator 
Tr(F^ u ) was already considered in [45], where it was found, with this ansatz, 
the following beta function: 



We would like to give an explanation of this result in the light of our holo- 
graphic effective action. In the case that the Wilsonian effective action was 
truncated to the lowest order operator, we should compare the result in [45] 
with ours without including the higher order operators. In this case our 
holographic effective action is only and exactly one loop as in the super- 
symmetric case. At first sight we have a discrepancy, since our result is only 
one loop while [45] exhibits a NSVZ-like beta function [16]. However, fol- 
lowing the analogy with the super-symmetric case, we should remind that 
in [45] it has been computed the canonically normalised Wilsonian effective 
action, while we have computed the effective action with the Wilsonian non- 
canonical normalisation. We can pass from one to the other one by means 
of a Jacobian that takes into account the rescaling anomaly, following [46]. 
This Jacobian has been computed in the M = 1 super-symmetric case and in 
the super-symmetric theory it accounts for the link between the holomorphic 
beta function, that is only one loop, and the beta function for the canonical 
coupling. Oddly the analogous computation does not seem to have been ever 
performed for pure QCD. We fill here this gap. In fact our computation 
mimics the one in the super-symmetric case, that is therefore recalled below. 
In the gauge theory with M = 1 super-symmetry, without matter multiplets, 
logJ, the logarithm of the Jacobian for the rescaling anomaly, receives a 
contribution from the gluons, V, from the gluoinos, ip, -0, form the auxiliary 
scalar field, D, and from the scalar field that implements the gauge-fixing 
condition in the Feynman gauge, c: 




(75) 



log J = logg lim (Tr v (exp(- — (-A A 5 fiU + i2ad F )) 
-Trv(exp(-— 7/1 L> M )) - Tr^(exp(- — (7^)) 
+Tr D (exp(-^-(-A A )) - Tr c (exp(-^-(-A A ))) 



(76) 
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Quite obviously in the pure gauge case we get: 

logj = logg lim (Tr v (exp(-^—(-A A S^ + i2ad F )) 



-Tr c (exp(-^(-A A ))) (77) 



We have computed that in this case: 

log J = logg^j 1 - £ J Tr(F^) 

Now, since the canonical coupling, g c , and the non-canonical Wilsonian cou- 
pling, gw, are related by the relation: 

hr "k" Pjl099 < m 

and the Wilsonian coupling is only one loop within our accuracy, we get for 
the canonical beta function, within our accuracy: 

^--^od-^A)- 1 (80) 

with: 

a 1 11 



167T 2 3 

A = T^l (81) 

in perfect agreement with the result in [45]. From Eq.(80)-(81) we can read 
that the contribution of the rescaling anomaly to fli in the pure gauge theory 
is: 

to = M> = (TeW? (82) 

where our convention is that positive f3o,Pi give origin to a negative beta 
function. Incidentally this seems to answer a question raised by Shifman [16] 
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about the operator versus canonical conformal anomaly in the pure gauge 
theory (within our accuracy, since in principle there could be contributions 
of higher order operators also in the Jacobian J) . This shows that the missing 
term that must come from higher order operators to get agreement with the 
two- loop perturbative result: 

Pl = Jvby f (83) 

is: 

= "(leM < 84 > 

We leave for the future the evaluation of higher order operators in our holo- 
graphic effective action. 



7 Miscellanea 



In this section we collect some comments about analogies with the existing 
literature. In recent years there has been much progress in understanding 
the non-perturbative physics of super-symmetric gauge theories based on 
such concepts as effective action, holomorphy, holography and integrable 
systems. One of the first non-trivial results was the exact NSVZ (canonical) 
beta function [16] in M = 1 super-symmetric theories, that is related to 
the fact that the Wilsonian (holomorphic) coupling constant gets only one- 
loop divergences because of the holomorphic properties of the Af = 1 SU SY 
action written in terms of super-fields. As we have just seen this distinction 
plays a role in this paper, if we want to relate our Wilsonian holographic 
effective action to the perturbative canonical one by a computation analogous 
to the one in [46]. Some time ago a holographic correspondence was found 
in [10], between a boundary M = 4 four dimensional super-symmetric gauge 
theory and a bulk five dimensional super-gravity (super-string) theory. New 
techniques based on integrable systems and holomorphic matrix models have 
been developed to find the low energy effective action of M = 2 and M = 1 
super-symmetric gauge theories in [47] and [48] respectively. The extension 
of these techniques to non-supersymmetric gauge theories in four dimensions 
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such as QCD4 encounters considerable difficulties. Though at technical level 
the construction in this paper of the holographic effective action in large- N 
QCD4 for planar self- avoiding loops seems peculiar, at least for the moment, 
to the M = theory, we would like to interpret the existence of such effective 
action in the light of some of the concepts that played such an important role 
in the super-symmetric case. We have already stressed that the construction 
of the effective action from the loop equation is a case of holography, i.e. of 
a correspondence between a boundary theory supported on loops, defined by 
the loop equation, and a bulk theory defined by the effective action. The 
key point for the existence of this holographic correspondence is holomorphy, 
that enters here in changing variables in the loop equation in such a way 
that the quantum contribution to the loop equation in these new variables 
involves the contour integral of the Cauchy kernel and thus can be reduced to 
the computation of a regularized residue. In this respect the Cauchy theorem 
can be considered the oldest and simplest case of holography. Yet, there are 
two other points of contact with the holography in the sense of [10]. The first 
one is that to make the quantum residue vanishing is necessary to map the 
local bulk degrees of freedom of the theory into the boundary cusps. Thus we 
have a correspondence between the degrees of freedom in the bulk and at the 
boundary, that is the distinguishing feature of holography in modern sense. 
In addition this is achieved by means of a conformal mapping that induces a 
new metric on the fundamental domain in the upper half plane, in analogy 
with the AdS theory [10]. Notice also that in our approach, in the four 
dimensional theory, there are logarithmic corrections to such metric due to 
the conformal anomaly, as already suggested in [35]. It is perhaps remarkable 
that these features arise directly from first principles, i.e. directly from the 
loop equation. 

The usual unitary monodromy of the connection that enters the loop equation 
is embedded into a non-unitary family of monodromies whose curvature is 
either of ASD or of SD type: this non-unitary family plays the analogous 
role of the holomorphic chiral ring in M = 1 super-symmetric gauge theories 
[49]. In fact the structure of the "classical" holographic action for large- ./V" 
QCD^ itself has many analogies with the one of N = 1 super-symmetric 
gauge theories [48]. This effective action is the sum of three terms. The 
classical action, the Veneziano-Yankielowicz [50] potential (essentially the 
one-loop contribution) and the contribution of the Konishi anomaly, that 
occurs as an anomaly in the conformal rescaling of the chiral super-fields 
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in the functional integral. It is interesting to observe that our "classical" 
holographic action has in fact a very similar structure. The Af = 1 effective 
action is defined as a functional of the glueball super-field, that is a scalar 
composite operator. In our Af = case the "classical" holographic action is a 
functional of a composite operator as well, the self-dual and anti-selfdual (in 
the Z 2 orbifold version) components of the curvature tensor, that play the role 
of the chiral curvature super-fields W and W in the super-symmetric case. 
There is even a hyper-Kahler structure, often associated to super-symmetry, 
that is inherited by the moduli fields of the Hitchin systems determined by 
these ASD or SD fields. Also the "classical" holographic action consists of 
three terms: the classical action, the localisation determinant (that is only 
one loop) and the Jacobian of the change of variables to the holomorphic 
gauge. This third term is the analogue of the Konishi anomaly and it is 
indeed a Jacobian under field-dependent gauge transformations living in the 
complexification of the gauge group. The occurrence of Hitchin systems in the 
holographic effective action can be considered the analogue of the occurrence 
of Hitchin systems in the solution for the low energy effective action of the 
Coulomb branch of Af = 2 super-symmetric gauge theories [47] . We would 
like also to recall the reader that the idea of using Hitchin systems to get 
control over the large- N limit of QCD arises by an energy-entropy argument 
in [51], that is rooted in the idea of dominating the large- N limit by the 
saddle-point method. Hitchin systems, being completely integrable, admit 
an " Abelianization" map, not to be confused with the one of section 3, that 
maps parabolic SU(N) Hitchin bundles into U(l) bundles over branched N 
coverings. Thus the system gets Abelianized in such a way that the number of 
local moduli passes from order of iV 2 for the parabolic Hitchin system to order 
of N for the Abelian system onaiV covering. Hence for large- N the entropy 
of the functional integration is re-absorbed into the Jacobian of the change 
of variables of the Abelianization map and the usual saddle-point argument 
for vector-like models applies. In fact we could interpret the holographic 
map as an attempt to understand in a more precise way, i.e. directly from 
the loop equation, how Hitchin systems help to get control over the large- TV 
limit of QCD. Finally, the use that we have made of the conformal map 
to a cuspidal fundamental domain resembles geometric engineering [52] and 
in particular the fact that only the local behaviour close to the singularity 
seems to count. The use of cuspidal singularities resembles also the idea [35] 
of considering five dimensional strings with the four dimensional boundary 
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at a singular point of the metric in the fifth coordinate, as required by the 
zig-zag symmetry. In turn the zig-zag symmetry plays a crucial role in the 
vanishing of the quantum residue in the loop equation. 



8 Conclusions 



We here advocate a drastic change in our way of attempting a solution of 
the loop equation in the large- N limit of QCD. Rather than trying to solve 
this equation on the entire loop algebra, a highly non-commutative problem, 
perhaps of non-hyperfinite nature, we restrict ourselves to the commutative 
hyper-finite algebra generated by a fixed self-avoiding loop. We look at this 
problem as a problem in holography, following the analogy with the Cauchy 
theorem, that reconstructs from the values of a holomorphic function on a 
loop its values in the internal region. We have constructed a holographic 
quantum effective action for the eigenvalues of the curvature of a twisted lo- 
cal system, i.e. a centrally extended Hitchin system. Hitchin systems occur 
in the functional integral by means of the localisation on a microcanonical 
ensemble labelled by the levels of the ASD (SD) curvature of the gauge 
connection. We have obtained a holographic form of the loop equation by a 
change of variable from the connection to the ASD (SD) curvature in the 
holomorphic gauge. We have written the holographic loop equation for a 
special connection, whose curvature is of ASD (SD) type. In the new in- 
tegration variable, i.e. the ASD (SD) curvature in the holomorphic gauge, 
the quantum contribution to the loop equation is reduced to the compu- 
tation of a regularized residue, that turns out to be loop independent for 
self-avoiding loops. To this form of the loop equation we have associated a 
holographic "classical" effective action. Finally, we have obtained the holo- 
graphic quantum effective action from the "classical" one by means of a 
conformal mapping of the loop equation to a cuspidal fundamental domain 
on which the regularized residue for self-avoiding loops vanishes identically. 
The conformal map makes possible a peculiar gauge choice and a consequent 
reduction of the curvature to a Borel sub-algebra. The physically interesting 
unitary Wilson loops are embedded into the algebra generated by Wilson 
loops of ASD (SD) type via a limiting procedure that involves a twisting of 
the Higgs field of the ASD (SD) Hitchin system. 
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In this paper we have computed the first coefficient of the beta function that 
arises from the holographic "classical" effective action for the connection in 
the unitary limit, finding exact agreement with the perturbative one-loop re- 
sult. In addition we have computed that part of the second coefficient which, 
following the analogy with the rescaling anomaly in M = 1 super-symmetric 
gauge theory, arises from the rescaling in the holographic effective action, in 
passing from the Wilsonian to the canonical coupling constant. We found 
also in this case exact agreement with an independent result based on con- 
tinuum exact renormalization group. The complete second coefficient of the 
beta function of the holographic effective action is left as a problem for the 
future. 
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